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I:  INTRODUCTION: 

Consider  an  insect  that  displays  observable  stages  throughout  its 
lifetime.   A  common  problem  is  to  estimate  the  average  or  typical  time 
required  for  the  insect  to  reach  a  particular  stage  (e.g.,  the  mature 
or  adult  stage)  of  development.   One  approach  to  this  problem  is  to 
observe  a  cohort  of  such  organisms  periodically  over  time  and  record 
the  number  of  insects  in  the  various  stages  and  average  the  times 
obtained. 

It  is  not  always  possible  to  observe  a  population  over  time 
since,  in  some  cases,  the  insects  can  not  be  recaptured  to  ascertain 
the  difference  between  previous  and  current  observing  time.   In  other 
cases  the  insects  may  have  to  be  killed  or  their  habitat  may  have  to 
be  destroyed  to  determine  which  stage  they  are  in. 

We  thus  take  samples  from  the  population  of  insects  at  a  sequence 
of  times.   We  let  t^,  t^ t^  denote  the  times  at  which  we  select 

a  sample.   The  samples  will  be  taken  randomly  and  without  replacement, 
although  we  will  assume  that  the  population  is  large  enough  that  the 
provision  for  sampling  without  replacement  will  not  affect 
calculations  and  the  samples  then  will  be  assumed  to  be  independent. 

Boyer  and  Deaton  (1984)  and  Pontius  (1987)  derived  and  evaluated 
nonparametric  estimators  for  the  average  time  of  development  across 
various  stages,  under  the  assumption  of  no  mortality  in  the  course  of 
the  experiment.   Here  we  extend  their  results  to  the  situation  with 
only  one  stage,  i.e.,  from  non-maturity  to  maturity,  but  where  the 
sampled  insects  might  be  either  dead  or  alive.   In  this  case,  death 
could  occur  either  before  or  after  maturity.   We  take  independent 
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samples  of  fixed  size  and  record  the  numbers  of  both  mature  and  non- 
mature  organisms,  noting  mortality  in  each  group,  at  various  observing 
times  and  estimate  the  average  time  to  maturity,  and  its 
corresponding  standard  error.   We  compare  the  mean  with  its 
theoretical  value  and  make  recommendations  for  appropriate  estimation 
procedures . 


II:  THE  RELATIONSHIP  OF  OBSERVATIONS  TO  THEORETICAL  QUANTITIES: 

Consider  an  experiment  in  which  we  are  interested  in  progression 
of  an  organism  from  one  stage  to  another  stage,  i.e.,  from  the  stage 
we  label  "  non-mature"  to  the  stage  we  label  "mature".   In  the 
performance  of  this  experiment,  we  collect  independent  samples  of 
insects  of  size  n^  at  time  t^,  i=l,2 F.   The  insects  collected 

may  be  either  dead  or  alive,  mature  or  non-mature.   We  then  allocate 
the  elements  of  the  ith  sample  to  one  of  four  categories  which  are 
non-mature  and  dead  (ND) ,  non-mature  and  alive  (NA) ,  mature  and  dead 
(MD) ,  mature  and  alive  (MA).   In  each  category  at  time  i,  the 
observations  have  the  binomial  (n.  ,  p..  )  distribution,  j-=  ND,  NA, 


MD,  MA,  and  the  overall  data  have  the  multinomial  (n     d      d 


IJ 


•  iND 


iNA  ' 


iMD 


p.wA^)  distribution. 


We  let  f(x,y)  represent  the  bivariate  density  of  X  (time  to 
maturity)  and  Y  (time  to  death).   We  presume  that  f(x,y)>0  only  if 
both  X  and  y  are  positive.   From  the  figure  2.1,  we  can  see  that  at 
time  t^  the  probabilities  for  the  categories  are  represented  as 

outlined  below. 


ND=  the  area  of  FBOX 


NA-  the  area  of  CBF 


MD-  the  area  of  Bt.  0 


iND 


iNA 


iMD 


t. 

ly 


t. 


t. 


f(x,y)  dx  dy. 


f(x,y)  dx  dy, 


iy*  ix 


'iy 


f(x,y)  dx  dy. 


MA=  the  area  of  Yt.  BC, 

ly 


iMA 


IX 


f(x,y)  dx  dy. 


Y  -  time  to  death 


IX 


X  -  time  to  maturity 


t.  -  the  shorthand  notation  for 


point  (t^,  0) 


t.  -  the  shorthand  notation  for 


point  (0,  z^ 


t.  -  time  of  ith  sample 

FIGURE  2.1 
Note  that  n^  is  the  size  of  the  sample  taken  at  the  ith  sampling 


time.   We  let  p.  -  p.„.  +  p. 


iMD 


t. 

IX 


f(x,y)dy  dx 


Consider  an  insect  that  displays  observable  stages  throughout  its 
lifetime.   Then,  again  from  figure  2.1,  we  let 

A 

p^  -  the  observed  proportion  of  insects  at  time  i  which  have 
reached  the  mature  stage  at  time  i 


#MA.+  #MD. 
1     1 


n 

1 


Here  we  let  #MA  be  the  number,  from  the  sample  of  size  n   found 
i.  i' 

in  the  MA  category,  and  #MD^  represents  the  corresponding  number  in 


the  MD  category  at  time  t. . 
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III:  THE  ESTIMATION  PROCEDURE  AND  THE  THEORETICAL  CALCULATION: 
a:  Expected  time  to  maturity 

Let  X  represent  time  to  maturity  and  let  Y  represent  time  to 
death.   Since  we  are  primarily  interested  in  the  length  of  time  which 
it  takes  to  reach  maturity,  we  restrict  our  attention  to  those 
organisms  that  live  long  enough  to  do  that.   Hence  the  parameter  that 
we  shall  estimate  is  the  average  time  to  maturity,  among  those  that  do 
reach  maturity.   Formally,  the  parameter  of  interest  is  T  =  E[X|X<Y] . 

Then,  one  can  show  that 

f(x|x<y)  -  f(x|x-y<0) 

-_f_  P(X<x|X-Y<0) 

ax 

a   P(X<x  and  X-Y<0) 


ax     P(X-Y<0) 

a  J-"  J  "  f(s,t)dtds 


s 


ax 


^  f(s,t)dtds 


t)dt 


If 


g  f(s,t)dsdt 


Hence 


T  -  E(X|X<Y)  -  j  I   xf(x|x<y)dx 

I   00   [   00 

J  0  J  x  ^f(x,y)dydx 


lol 


00 


X 


f(x,y)dydx 


equivalently,  using  change  of  order  of  integration, 


_  J  0  J  0  xf(x,y)dxdy 

[  "  I  y  f(x,y)dxdy 

J  0  J  0  (3.a.l) 

Note  that  we  assume  the  density  to  be  zero  when  either  argument 
is  negative,  otherwise  the  lower  limits  in  the  expression  above  would 
all  be   -  00. 

b:  An  estimator  of  T 

In  the  previous  expression  for  T,  let  us  denote  the  numerator  by 
U  and  the  denominator  by  V.   Then  note  that 

^  "  J 0  J 0  xf(x,y)dxdy  =  J^  J^  xf(x.y)dydx. 
If  the  X-axis  is  broken  into  inteirvals 

[O.t-j^]  ,  [t^,t2]  ,  [t2,  t^]  , then  clearly  U  is  reasonably  approximated 

(in  a  manner  similar  to  Riemann  Sums)  by 

t,  A    A       t,+  t,  A    A       t-+  t„ 

U  =  Pi  ( )  +  (  P2-Pi)(_J i_)  +  (P.-Po)(_J 1_)  + 

2  '"2  ^     ^  2 

A        A  t„+   t„   ^ 


^F  ^'F-1' 


(See  Figure  3.1  below,  and  recognize  that  p   is  the  fraction  of  the 

sample  which  lies  above  the  line  X  -  Y  and  to  the  left  of  t   ) 

ix'  ■ 

In  a  similar  manner  V  is  just  the  fraction  of  the  population 

A 

which  lies  above  the  line  X  =  Y,  and  is  approximated  by  p„,  the 

F 

fraction  of  mature  organisms  found  at  the  last  sampling  time.   Of 

course,  this  presumes  that  sampling  time  F  is  sufficiently  large  that 

nearly  all  organisms  that  will  mature  have  done  so. 
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AAA 


Thus  we  arrive  at  an  estimator  T  by  taking  T  =  U  /  V,  this  gives 


us 


*  *        *^1  A        ^  t„+   t^  A        A  t-+   t_ 

T  =  [  Pl( )  +  ( p,-  Pi)  (_ _)  +  ( p,-  p,)  (^ i 


'2  T 


'3  ^2^ 


)  + 


•••■-^  (  P..-P..o)(_Zi_J_l_)+  (  Pp-  Pp.i)Pp    ^   ^-^  )]/Pp 


'F-  ^F-2' 


(3.b.l) 
and  with  a  little  simplification  we  get  an  equivalent  form  which  is 
slightly  better  for  some  computations. 


^    Pi      F-1   t..^  -  t.^^     p 


'2    ^1 

_  (     )  +  2    ( 


i-2 


)  (^)  ,  (  "^  ^   'f-1  ^ 


Y  =  time  to  death 


1        i+1 
X  -  time  to  maturity 

t.-  sample  time  ,  i-1,2 F 

FIGURE  3.1 
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A 

c:  Approximate  mean  of  T 

Using  the  approximate  formula  (  Mood,  Graybill  and  Boes,  1974), 
we  find 

E(_^)«  _J_   -        ^        Cov(A,B)  +  _J_   Var(B). 
B     Mk       ^  3 

Hence  we  note  that  after  a  little  simplification  (because 

A      A 

cov(p^,Pp)  =  0  ) 

A 

,  Pi  .   Pi  ,  "fPf^i-Pf  ,   ^  Pi 

^'^—^ >   ~  ( )  --->  as  n   --->co  . 

Pp      PP      "pPp  Pp 

Pi 
Therefore   ^   is  an  approximately  unbiased  estimator  of  the 

Pp 
corresponding  ratio  of  population  parameters  as  n_--->  «,  and,  since  T 
is  nearly  a  weighted  average  of  such  quantities  we  are  led  to  believe 

A 

that  T  will  be  unbiased  or  nearly  unbiased  as  an  estimator  of  T. 

A 

d:  Approximate  variance  of  T 

Again  we  use  an  approximation  (Mood,  Graybill,  and  Boes,  1974) 

Var(T)  -  Var(  ^  ) 
B 

2 


/^ 


A  ,  Var(A)    Var(B)      2Cov(A,B)  ,  ,.   ^   ^^ 
-s— ^ —  ^  * ^J  (o.d.l). 


'b     ^k  ^Y>  '^A^B 


In  this  case 


-  ^2  "  ^-2   t   -  t     .     t„+  t„  .  . 

^  -  "  -  -A  Pi  ^  -    (     '    '^'  )Pi  ^(   ^   "-')Pp. 
2         i-2        2       ^       2      ^ 


A        A 
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B  =  V  =  Pp. 

to     F-1   t,  .  -  t, ,,       r+  t 


Thus  ,^  =  [  C_2)V,  -,l^(.±}_J^)V,H_L_^)v/, 


^B  -  Pf  • 


-  ^2   2   Pid-Pi)   F-1    t    -  t.  ^   2  p.(l-p.) 
Var(A)  =  ( f_)     ^    ^   +2    (   ^'^    ^^^  )   ^   ^ 

2         n        i-2       2  n. 

J-  1 

^(  V  ^F-1  ^2  Ppd-Pp) 


■^F 


Var(B)  =    ^    ^ 


"f 


Cov(A,B)  =  (   ^   ^'^)  Var(pp) 


_  (  V  ^f-1  ^  Pf<^-Pf> 


np 


In  estimating  the  variance  of  T,  we  will  insert  the  corresponding 

A 

sample  quantities  (p!^s)  into  each  term  and  replace  in  (3.d.l). 
e:  The  comparison  to  the  estimator  of  Pontius 

A 

The  p^  we  have  defined  represents  the  proportion  of  organisms 

A 

sampled  at  time  i  that  have  reached  the  mature  stage,  so  that  1-p. 

A 

(i.e.,q^)  represents  the  proportion  that  have  not  yet  reached  the 
mature  stage. 
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^  ^1      ^  F  A         A  t.+   t.    ^ 

Since  T  =  [ p  +  2   (  p  -  p  ,)(   ^   ^'^  )1 

2   ^   i=2    ^   ^'^  2 


Pp 


and  p .  =  1 - q -  . 

tl         F-1      A  t.    ,-    t.    ,  A         t   +   t 


2    i=2 


T  = 


i-qi 


Under  the  assumptions  made  by  Pontius,  t^  must  be  zero  and  p 

U  F,  S 

A 

A 

(our  qp)  must  be  zero.   Using  this  condition  the  expression  for  T 
reduces  to 

t      F-1  A    t.  ,-  t.  , 
T-^+  S   q.(   ^-^^    ^-^  )■ 
2     i-2  ^     2 

This  is  exactly  the  same  as  the  estimator  derived  by  Pontius 

(1987)  for  the  case  in  which  it  is  assumed  that  no  insects  have  died. 

Hence,  although  we  have  extended  the  method  for  computations  to 

populations  which  may  contain  dead  insects,  the  final  result  is  still 

the  same  as  the  result  derived  by  Pontius  if  no  mortality  is  observed 

during  the  experiment. 
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IV:  NUMERICAL  INTEGRATION  AND  THE  SCOPE  OF  THE  SIMULATION  STUDY: 
a:  The  numerical  integration 

Since  the  quantity  T  =  E(X|X<Y)  (which  is  the  parameter  of 
interest  in  this  problem)  is  not  readily  derived  for  the  bivariate 
normal  distribution,  we  used  the  IMSL  subroutine  DBLIN  to  numerically 
perform  the  two  double  integrations  required  (expression  3.a.l).   A 
copy  of  the  numerical  integration  program  is  in  the  appendix, 
b:  The  scope  of  the  simulation  study 

We  developed  a  computer  program  to  simulate  the  experiment 
described  above.   As  an  underlying  distribution  we  used  the  bivariate 
normal  distribution,  which  has  a  large  number  of  shapes  depending  on 
one's  choices  of  parameters.   Since  we  wanted  to  investigate  the 
lifetimes  of  insects  which  mature  before  death,  and,  since  in  most 
populations,  lifetime  is  much  longer  than  maturation  time,  we  chose 
mean  of  "time  to  death"  (i.e.  u)    to  be  larger  than  mean  of  "time  to 

maturity"  (i.e.  n^   ).   For  normal  distributions,  a  distance  of  three 

standard  deviations  from  the  mean  in  both  directions  will  cover  most 

of  the  observations.   Hence  we  arbitrarily  chose  u     and  u   and  then 

X     '^y ' 

picked  standard  deviations  of  "time  to  death"  and  "time  to  maturity" 

(i.e.  a       a       resp.)  to  be  values  such  that  u  ±  3a  and/or  u   -    3a 
J         '^  X    X    '    '^y    y 

were  in  between  the  first  and  the  last  observation  time.   (We  were 

willing  to  allow  some  of  the  insects  to  die  after  our  experiment  was 

over;  that's  why  we  are  not  concerned  with  n  +3a     being  inside  our 

observation  time.)   We  let  n     be  fixed  at  6  and  a     be  fixed  at  2 .   We 

X  X 

then  let  fi^-  9   with  a^-   2  and  a^-   3,  u-   12  with  a  -  3  and  a  -  4  and 


15 
y  =  15  with  a  =  4  and  a  =■  5.   We  also  allowed  the  correlation 

coefficient  p  to  be  0,  0.25  and  0.5.   We  considered  samples  of  size  15 
and  30  at  each  time  period  and  we  had  sampling  regimens  which  had  10 
sampling  times  equally  spaced  (every  two  time  units  from  time  t  =  2  to 
time  t  =  20)  and  unequally  spaced  (also  from  t  =  2  to  t  =  20) .   For 
each  of  the  72  parameter  configurations,  100  repetitions  of  the 
experiment  were  run,  and  computations  were  made. 

A 

For  each  repetition,  T  and  its  estimated  standard  error,  which 
will  be  denoted  by  std(T) ,  were  calculated.   Additionally  the 

^  A  A  A 

intervals  T  ±  2std(T)  and  T  ±  3std(T)  were  computed. 

For  each  parameter  configuration,  several  summary  statistics  over 
the  100  repetitions  were  retained  and  examined.   They  include: 

A 

(1)  Average  value  of  T. 

A  A 

(2)  Number  of  intervals  of  the  form  T  ±  2std(T)  that 
contained  T. 

A  A 

(3)  Number  of  intervals  of  the  form  T  ±  3std(T)  that 
contained  T. 

A 

(4)  Sample  standard  deviation  of  T. 

A 

(5)  The  interval  for  the  mean  of  T,  given  by  average  value 

*  A 

of  T  ±  2  *  sample  standard  deviation  of  T. 

A 

(6)  Average  value  of  estimated  standard  error  of  T. 

A  copy  of  the  simulation  program  is  found  in  the  appendix, 
c:  Two  possible  computations: 
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The  population  p^  in  formula  (3.b.l)  should  be  monotonically 

increasing  as  a  function  of  i,  but  when  we  dealt  with  the  sample 

proportions,  we  occasionally  found  p.  to  be  smaller  than  p.   .   We 

used  both  the  negative  value  and  zero  respectively   for  the  quantity 

P^-  P^.i  i"  °^^  simulation  study  when  p.  was  smaller  than  p.  -.  . 
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V:    CONCLUSIONS: 


A        A 


We  discovered  early  that  the  quantity  p.-  p.   would  be  negative 
sometimes,  even  though  the  corresponding  population  quantity  is  not. 

A        A 

We  considered  replacing  these  negative  values  of  p.-  p.   by  zero  in 

the  computations.   In  fact  in  early  stages  of  the  simulation,  we  did 
computations  both  ways.   After  comparing  the  results,  we  recommend  the 
use  of  the  negative  value  instead  of  using  zero  when  this  occurs.   Our 
recommendation  is  based  on  the  comparison  of  the  simulation  results  to 

A 

the  numerical  integration  results,  i.e.,  the  average  T  from  100 

A        A 

simulations  with  p.-  P^_-,  allowed  to  be  negative  is  considerably 

A        A 

closer  to  the  numerical  intergration  than  p.-  p.    set  to  zero.   That 
is,  the  average  T  from  100  simulations  tends  to  be  biased  with  p  - 

A 

p.  ,  set  to  zero. 

Using  the  computer  output  of  the  simulation  program,  we  have 

A 

obtained  an  estimated  value  of  T  for  each  different  configuration  of 
Mj^.  M  ,  CT^,  CT  ,  p.   We  have  calculated  T  as  the  average  of  the 


K  A 


estimated  T's  and  formed  an  interval  of  the  form  T  ±  Iok    .   Since  we 

T 

A 

have  100  independent  T's  from  100  independent  repetitions,  by  the 

Central  Limit  Theorem  T  will  be  apporximately  normally  distributed  and 

therefore,  the  quantity  T  ±  2a ^    is  approximately  a  95%  confidence 

T 
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A 

interval  for  E(T) .   We  found  that  for  (1)  29  out  of  36  parameter 
configurations  with  sample  size  15.  (2)  26  out  of  36  with  sample  size 
30,  (3)  27  out  of  36  with  equally  spaced  sampling  time,  (4)  28  out  of 
36  with  unequally  spaced  sampling  times, the  confidence  intervals 
contained  the  parameter  of  interest  (T  or  E[X|X<Y]). 

A 

Note  that  the  standard  error  of  T  is  estimated  in  two  ways  in 
this  study.   The  first  is  from  looking  at  100  independent  repetitions 
of  the  experiment  and  estimating  the  standard  deviation  of  the 

A 

population  of  T's.   This  quantity  is  used  in  the  confidence  interval 
calculations  described  above.  The  second  method  is  from  the  series 
expansion  for  an  individual  repetition,  and  uses  formula  (3.d.l).   We 
have  averaged  these  values  and  compared  them  with  the  lengths  of  the 
confidence  intervals  described  above.   (Note  that  the  length  of  the 

A 

confidence  interval  should  be  4  x  a^  /  J  100  or  approximately  0.4 

A 

times  the  standard  error  of  T) .   We  find  good  agreement  with  this 
criterion  in  our  tables,  indicating  that  the  short  series  expansion 
gives  a  good  estimate  of  the  actual  standard  error. 

Note  also  that  the  standard  errors  for  samples  of  size  15  are 

approximately  1.4  (or  J2  )  times  the  corresponding  standard  errors  for 
samples  of  size  30.   This  is  an  indication  to  us  that  even  though 
there  is  no  such  term  directly  found  in  the  Taylor  series  expansion 
which  gave  us  (3.d.l),  (at  least  for  the  situation  where  all  the 
samples  are  of  the  same  fixed  sample  size  n) ,  the  standard  error  of 

the  estimator  may  be  proportional  to  1/Jn  ,  as  are  standard  errors  of 
means . 
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Another  study  that  we  performed  is  summarized  in  table  2.   There 
we  used  the  standard  errors  for  each  repetition,  using  the  formula 

A 

(3.d.l)  and  denoted  here  by  "std(T)"  ,  and  formed  the  intervals 

A  A  A  A 

T±2std(T)  and  T±3std(T) .   Then  we  determined  how  many  of  the  100 
repetitions  generated  confidence  intervals  which  contained  the 
parameter  values.   Although  the  results  do  not  show  perfect  agreement 
with  the  nominal  coverage  proportions,  they  always  show  around  80  to 
95  percent  of  confidence  intervals  which  contained  the  parameter 
values . 

A 

As  a  nonparametric  estimator,  it  was  not  possible  to  show  that  T 
is  a  Uniform  Minimum  Variance  Unbiased  Estimator  (UMVUE) ,  but  we  did 
demonstrate  that  it  is  nearly  unbiased  for  the  parameters  and 
distribution  used  in  this  study.   We  also  showed  that  confidence 
interval  coverage  is  near  the  nominal  levels. 

From  the  simulations  we  have,  we  did  not  find  any  significant 
difference  in  results  between  10  equally  spaced  and  unequally  spaced 
observation  times.   It  seems  to  be  a  convenience  for  an  experimenter 
not  to  have  to  follow  a  strict  schdule  to  observe  and  collect  data. 

The  preceding  work  is  not  a  complete  study  of  the  proposed 
estimator.   Future  work  could  fill  several  gaps.   Foremost  among  these 
might  be  a  study  using  alternative  underlying  bivariate  distributions. 
(The  bivariate  normal  was  used  in  this  case  because  of  ready  access  to 
it,  the  ease  of  computations,  and  the  versatility  in  terms  of  the 
number  of  shapes  availables.)  Other  aspects  which  should  be 
investigated  include  a  more  thorough  study  of  the  effect  of  sample 
size,  a  study  of  the  effects  of  unequal  sample  sizes  (at  the  various 
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observation  times),  and  a  more  thorough  study  of  the  lengths  of  the 
intervals  between  sampling  times  (how  big  should  they  be,  and  do  they 
need  to  be  equal  ?). 


u: 
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APPENDIX 

REAL*8  DSEED,UX,UY,VX,VY,PXY,A,B,C,D,DEXP,DSQRT,FD,PP,X,YCD,CN 
*, AREA, FN, AY, F, BY 

DIMENSION  SIGMA(3) ,XYMA(200, 2) ,WKVEC(2) ,KKVEC(2) ,T(30) ,MD(30) 

,MA(30) ,ND(30) ,NA(30) ,P(30) ,F(30) ,G(30) ,G0(30) ,ET(100) 
*  ,ET0(100) 

COMMON   UX,  UY,  VX,  VY.PXY 
INTEGER   lER 
EXTERNAL  FN,  FD,  AY,  BY 
C 

c      *************************************** 
c     *  * 

C      *    T(i)  -  SAMPLE  TIME  AT  TIME  i     * 

C      *       i  =  l,2 10  * 

C      *    AX,  BX  ARE  LOWER  AND  UPPER  LIMIT  * 
C      *      ON  X  * 

C      *  * 

c      *************************************** 

AX-O . DO 

BX-30 . DO 

AERR-0. 000001 
C 

c     ******************************************* 
c    *  * 

C     *      UNEQUALLY  SPACED  OBSERVATION  TIME  * 
C     *         ALTERNATIVELY  * 

C     *       SET  T(i)  -  2  TO  20  BY  2  IS        * 

C     *         EQUALLY  SPACED  OBSERVATION  TIME  * 
C     *  * 

c     ******************************************* 

c 

T(1)-2.D0 

T(2)-3.DO 

T(3)-6.D0 

T(4)-9.D0 

T(5)-10.D0 

T(6)-11.D0 

T(7)-12.D0 

T(8)-15.D0 

T(9)-16.D0 

T(10)-20.D0 

DSEED     -  466364003. DO 

NR        -  30 

K        -  2 

IR        -  200 


23 

KKVEC  (1)  =1.D0 

WKVEC  (1)  =  O.DO 
C 

c  ************************************** 

c  *  * 

C      *   READ  18  PARAMETER  CONFIGURATION  * 
C      *  * 

c  ************************************** 

c 

501  READ(5,105,END=500)  UX.UY.VX.VY.PXY 

105  FORMAT (5F5. 2) 
C 

C  ************************************************* 

c  *  * 

C      *     SET  THE  FOLLOWING  VARIABLES  TO  BE  ZERO    * 
C      *     NSAM  -  #  OF  SAMPLE  SIZE.  * 

A 

SUM  =  SUM  OF  T.  * 

A 

*  SUMV  =  SUM  OF  STD(T).  * 
"it                                                                                    ^ 

*  ETSQ  =  SUM  OF  SQUARE  OF  T.  * 


* 


* 


g  t 

8  * 
8 

C  *  V    =  CROSS  PRODUCT  OF  T  IN  ORDER  TO     * 

C  *  COMPUTE  THE  VARIANCE  OF  T.         * 

8  t  ETMIN  -  MIN  OF  T.                        * 

8  *  ETMAX  -  MAX  OF  T.                        * 

8  *  C0UNT2  -  COUNT  THE  #  OF  2STD(T)  INTERVALS  * 

C  *  THAT  CONTAIN  THE  PARAMETER  OF    * 

C  *  INTEREST.                       * 

8  t  COUNTS  -  COUNT  THE  #  OF  3STD(T)  INTERVALS  * 

C  *  THAT  CONTAIN  THE  PARAMETER  OF    * 

C  *  INTEREST.                      * 

C  *  W     -  CROSS  PRODUCT  OF  STD(T)  IN  ORDER  * 

*  TO  COMPUTE  THE  VARIANCE  OF  STE(T)* 

C  *  THE  ABOVE  VARIABLES  ARE  COMPUTED  WITH  THE   * 

f*  -Je  A      A 

g  I  QUANTITY  P^-  P^_^  SET  TO  NEGATIVE  , '  IN-    * 

C  *  STEAD  OF  THAT  SUMO,  SUMOV,  ETSQO,  ETOMIN  * 

C  *  ETOMAX,  C0UNTG2,  COUNT03 ,  WO,  VO  ARE     * 

§  I  COMPUTED  WITH  THE  QUANTITY  P^-  P    SET   * 

C  *  TO  ZERO.                               * 

C  * 

C  ************************************************* 


* 


* 
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NSAM=0 
SUM=0 
SUMO=0 
SUMV=0 
SUMOV=0 
ETSQ=0 
ETSQO=0 
V=0 
V0=0 

ETMIN=100 . DO 
ETMAX=0 . DO 
ET0MIN=100.D0 
ETOMAX=0 . DO 
COUNT2=0 . DO 
COUNT3=0 . DO 
COUN02=0 . DO 
COUN03=0.D0 
SUMV=0 . DO 
W=0 .  DO 
W0=0 .  DO 

CALL  ERRSET(208,  256,  -1,  1) 
C 

c      ****************************************** 
c     *  * 

C  *  DOUBLE  INTEGRATION  BY  USING  DBLIN  * 

C  *      SUBROUTINE.  * 

C  *  CN  -  NUMERATOR.  * 

C  *  CD  -  DENOMINATOR.  * 

C  *  AREA  =  E(X|X<Y).  * 

C  *  ^ 

c      ****************************************** 

c 

CN  =  DBLIN(FN,  AX,  BX,  AY,  BY,  AERR,  ERROR,  lER) 
CD  -  DBLIN(FD,  AX,  BX,  AY,  BY,  AERR,  ERROR,  lER) 
AREA  -  CN/CD 
WRITE(6,121) 


121 
* 


FORMAT('-M4X,'AY',9X,'BY',10X,'AERR',8X,'IER',11X,'CDM1X  'CN' 
,10X,'AREA') 


WRITE(6,122)  AX,  BX,  AERR,  lER,  CD,  CN,  AREA 

122  FORMAT('-M0X,7(F10.6,2X)) 
C 

C  ****************************************** 

C  *  ^ 

C  *      100  REPETITIONS.  * 

C  *      SET  TT  AND  VARA  -  0.  * 

C  *  * 

c  ****************************************** 

c 
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DO  6  KK=1,100 

TT=0 

TTO=0 

VARA=0 . DO 
C      *********************************************** 
C      *  * 

C  *  GERERATE  BIVARIATE  NORMAL  RANDOM  NUMBER  * 

C  *  NR  =  #  OF  OBSERVATIONS  * 

C  *  K  -  #  OF  RADAM  VARIABLE  * 

C  *  SIGMA  =  VARIANCE -COVARIANCE  MATRIX  * 

C  *               (UPPER  TRIANGAL)  * 

C  *  IR    =  #  OF  ROW  DIMENSION  * 

C  *  XYMA  =  OUTPUT  NR  BY  K  MATRIX,  * 

C  *  lER    =  ERROR  PARAMETER  (OUTPUT)  * 

C  *  * 

c  *********************************************** 

SIGMA(1)=VX 

SIGMA(2)=PXY 

SIGMA(3)-VY 

IF  (KK.GT.  1)  GO  TO  103 

CALL  GGNSM  (DSEED,  NR,  K,  SIGMA,  IR,  XYMA,  WKVEC,  lER) 
C 

c      ****************************************** 
c     *  * 

C  *  DETERMINE  THE  GERERATED  DATA  THAT  * 

C  *      BELONG  TO  WHICH  OF  THE  FOUR  * 

C  *      CATEGORIES.  * 

C  *  NA  =  NON- MATURE  AND  ALIVE.  * 

C  *  ND  -  NON-MATURE  AND  DEAD.  * 

C  *  MA  =  MATURE  AND  ALIVE.  * 

C  *  MD  -  MATURE  AND  DEAD .  * 

C  *  * 

c      ****************************************** 

c 

NA(1)-0.D0 
ND(1)-0.D0 
MA(1)-0.D0 
MD(1)-0.D0 
DO  1  1-1,15 
X  -  UX+  XYMA  (1,1) 
Y  -  UY+  XYMA  (1,2) 

IF  (X.LE.T(l)  .AND.  Y.LE.T(l))  GO  TO  5 
IF  (X.LE.T(l)  .AND.  Y.GT.T(l))  GO  TO  2 
IF  (X.GT.T(l)  .AND.  Y.LE.T(l))  GO  TO  3 
IF  (X.GT.T(l)  .AND.  Y.GT.T(l))  GO  TO  4 
5  IF  (X.GT.Y)  GO  TO  3 
MD(1)-MD(1)+1 
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GO  TO  1 

2  MA(1)=MA(1)+1 
GO  TO  1 

3  ND(1)=ND(1)+1 
GO  TO  1 

4  NA(1)-NA(1)+1 
1   CONTINUE 

K=l 

C 

C      ***************************************** 

c     *  * 

C      *    CALCULATE  P.,  F. ,  G.  AND  * 

1    1    1 

C      *  F  * 

C      *      TT=SF.*G.  * 

C      *  i-1   ^    ^  * 

A  A 

C      *      (T  -  TT/Pp  LABELED  BY  ET(KK))    * 

C      *  * 

c      ***************************************** 

c 

P(K)=(MA(K)+MD(K) )/15 . DO 

IF  ((K-l).EQ.O)  T(K-1)-0.D0 

F(K)=(T(K)+T(K-1))/2.D0 

IF((K-l).EQ.O)  P(K-1)-0.D0 

G(K)-P(K)-P(K-1) 

IF(G(K).LT.O)  THEN 

GO(K)-O.DO 

ELSE 

GO(K)-G(K) 

END  IF 

TT-TT+F(K)*G(K) 

TTO-TTO+F(K)*GO(K) 

GO  TO  109 
103  DO  7  LL-1,10 

L-LL 

GO  TO  108 
109   DO  11  L-2,10 
108   CALL  GGNSM  (DSEED,  NR.  K,  SIGMA,  IR,  XYMA,  KKVEC,  lER) 

NA(L)-O.DO 

ND(L)-O.DO 

MA(L)-O.DO 

MD(L)-O.DO 

DO  12  1-1,15 

X  -  UX  +  XYMA  (1,1) 

Y  -  UY  +  XYMA  (1,2) 

IF  (X.LE.T(L)  .AND.  Y.LE.T(L))  GO  TO  15 

IF  (X.LE.T(L)  .AND.  Y.GT.T(L))  GO  TO  25 
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IF  (X.GT.T(L)  .AND.  Y.LE.T(L))  GO  TO  35 

IF  (X.GT.T(L)  .AND.  Y.GT.T(L))  GO  TO  45 
15    IF  (X.GT.Y)  GO  TO  35 

MD(L)=MD(L)+1 

GO  TO  12 
25   MA(L)=MA(L)+1 

GO  TO  12 
35   ND(L)=ND(L)+1 

GO  TO  12 
45   NA(L)=NA(L)+1 
12   CONTINUE 

P(L)=(MA(L)+MD(L))/15.D0 

IF  ((L-l).EQ.O)  T(L-1)=0.D0 

F(L)-(T(L)+T(L-1))/2.D0 

IF  (  L.EQ.IO  )  GO  TO  16 

VARA=VARA+(((T(L-l)-T(L+l))/2)**2)*P(L)*(l-P(L))/15.DO 
GO  TO  17 

C 

c      ****************************************** 
c     *  ^ 

C      *    CALCULATE  VARA,  VARB,  COVAB,  UA,  UB  * 
C      *     BY  INSERTING  SAMPLE  QUANTITIES  IN  * 


ORDER  TO  CALCULATE  VAR(T) ,         * 


6 

C      *     LABELED  BY  VARET.  * 

C      *  * 

c  ****************************************** 

c 

16  VARA-VARA+(((T(L)+T(L-1))/2)**2)*P(L)*(1-P(L))/15.D0 
VARB-P(10)*(1-P(10))/15.D0 

COVAB-( (T(L)+T(L-1) )/2)*VARB 

17  IF((L-l).EQ.O)  P(L-1)-0.D0 
G(L)-P(L)-P(L-1) 

IF  (  G(L).LT.O)  THEN 
GO(L)=O.DO 
ELSE 

GO(L)-G(L) 
END  IF 

TT-TT+F(L)*G(L) 
TTO-TTO+F ( L) *G0 ( L) 
IF  (LL  .GE.  1)  GO  TO  7 
11  CONTINUE 
UA-(TT)**2 
U0A-(TT0)**2 
UB-(P(10))**2 
GO  TO  110 
7   CONTINUE 
UA-(TT)**2 
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U0A-(TT0)**2 

UB=(P(10))**2 

LL=0 
110  NSAM=NSAM+1 

ET(KK)  =TT/P(10) 

ETO(KK)  =TT0/P(10) 

VARET-(UA/UB)*( (VARA/UA)+(VARBAJB) - (COVAB/(UA*UB) ) ) 

VAROET=(UOA/UB)*( (VARA/UOA)+(VARB/UB) - (COVAB/(U0A*UB) ) ) 

STDET=SQRT(VARET) 

STDOET=SQRT(VAROET) 

SUMV-SUMV+STDET 

SUMOV=SUMOV+STDOET 
C 

c      ******************************************* 
c     *  * 

C  *     CALCULATE  2  &  3STD(T)  INTERVALS     * 

C  *      WITH  LABEL  DL  AND  DU  AS  LOWER      * 

C  *       AND  UPPER  LIMIT.  * 

C  *  * 

c  ******************************************* 

c 

DL2-ET(KK)-2*STDET 

DU2=ET(KK)+2*STDET 

DL3-ET(KK)-3*STDET 

DU3-ET(KK)+3*STDET 

D0L2-ET0(KK) -2*STD0ET 

D0U2-ET0 (KK) +2*STD0ET 

D0L3-ET0 (KK) - 3*STD0ET 

D0U3-ET0(KK)+3*STD0ET 

IF  (  AREA. LE.DU2. AND. AREA. GE.DL2)  THEN 

COUNT2-COUNT2+1 

ELSE 

COUNT2-COUNT2 

END  IF 

IF  (  AREA. LE.DU3. AND. AREA. GE.DL3)  THEN 

COUNT3-COUNT3+1 

ELSE 

COUNT3-COUNT3 

END  IF 

IF  (  AREA. LE.D0U2. AND. AREA. GE.D0L2)  THEN 

COUN02-COUN02+1 

ELSE 

COUN02-COUN02 

END  IF 

IF  (  AREA. LE.D0U3. AND. AREA. GE.D0L3)  THEN 

COUN03-COUN03+1 

ELSE 
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COUN03=COUN03 

END  IF 

SUM=SUM+ET(KK) 

SUMO=SUMO+ETO(KK) 

ETSQ=ETSQ+ET (KK) *ET (KK) 

ETSQO=ETSQO+ETO (KK)*ETO (KK) 

IF  (ET(KK).LT.ETMIN)  THEN 

ETMIN=ET(KK) 

ELSE 

ETMIN=ETMIN 

END  IF 

IF  (ET(KK).GT.ETMAX)  THEN 

ETMAX=ET(KK) 

ELSE 

ETMAX=ETMAX 

END  IF 

IF  (ETO(KK).LT.ETOMIN)  THEN 

ETOMIN=ETO(KK) 

ELSE 

ETOMIN-ETOMIN 

END  IF 

IF  (ETO(KK).GT.ETOMAX)  THEN 

ETOMAX=ETO(KK) 

ELSE 

ETOMAX-ETOMAX 

END  IF 

IF  (NSAM  .EQ.  1)  GO  TO  119 

V=V+(NSAM*ET(KK)-SUM)**2/(NSAM*NSAM-NSAM) 

V0-V0+(NSAM*ET0(KK)-SUM0)**2/(NSAM*NSAM-NSAM) 

W-W+((NSAM*STDET-SUMV)**2)/(NSAM*(NSAM-1)) 

W0-W0+((NSAM*STD0ET-SUM0V)**2)/(NSAM*(NSAM-1)) 
GO  TO  120 

119  V-O.DO 
VO-0 . DO 
W-O .  DO 
WO-O .  DO 

120  IF  (KK.LT.IOO)  GO  TO  6 
PER2-COUNT2/100 
PER3-COUNT3/100 
PER02-COUN02/100 
PER03-COUN03/100 
XBAR-SUM/100 
XOBAR-SUMO/100 
STBAR-SUMV/100 
STOBAR-SUMOV/100 
VAR-V/99 
VARO-VO/99 
VARST-W/99 
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VAR0ST=W0/99 
STD=SQRT(VAR/100) 
STD0=SQRT(VAR0/100) 
STDST=SQRT(VARST) 
STDOST=SQRT(VAROST) 
C 

c  ******************************************** 

c  *  * 

C      *     CALCULATE  THE  UPPER  AND  LOWER  LIMIT  * 

C         *  ;^    A  * 

C  *      OF  THE  INTERVAL  FORMED  T±2cr«  .      * 

C  *  T        * 

C  *  * 

C  ******************************************** 

c 

UPPER=XBAR+2*STD 
DL0W=XBAR-2*STD 
UPPER0=X0BAR+2*STD0 
DLOW0=X0BAR- 2*STD0 
WRITE(6,107) 

107      FORMATCIMOX,  'UX'  ,10X,  'UY'.IOX,  'VXMOX,  'VY'  ,  lOX,  'PXY' ) 
G 

c  ****************************************** 

c  *  * 

C  *  WRITE  THE  PARAMETER  CONFIGURATIONS* 
C  *  AND  THE  SUMMARY  STATISTICS.  * 
C      *  * 

c      ****************************************** 

c 

WRITE(6,106)  UX,UY,VX,VY,PXY 
106  FORMATC-' ,9X,5(F5.2,7X)) 
WRITE(6,100) 

100  FORMATC '-MOX.'NR',  5X, 'K' ,  2X,  '  SIGMAM8X,  '  IR' ) 
WRITE(6,101)  NR,  K, (SIGMA(J) , J-1 , 3) ,  IR 

101  FORMATC-',  8X,  14,  2X,  14,  2X,  3(F5.3.2X),  14/) 
WRITE(6,115) 

115   FORMATC '-',5X, 'THE  FOLLOWING  VALUE  IS  WHEN  THE  PROBABILITY  IS  NEG 
*ATIVE:') 

WRITE C 6, 104)  TT,  ETCKK) ,  SUM.XBAR, ETSQ.VAR, STD.ETMIN, ETMAX.DLOW 
*  UPPER , C0UNT2 , PER2 , C0UNT3 , PER3 , STBAR , VARST , STDST 
104  FORMAT  C  '  -  MOX ,  '  TT- '  ,  FIO .  5  ,  lOX ,  '  EXPECTATION  OF  TIME  TO  THE  MATURI 
*TY  -',  F10.5/10X,'SUM-',F10.5,10X,'XBAR-',F10.5  /lOX.'SUM  OF  SQUAR 
*E  -', FIO. 5, /lOX,' VARIANCE  -' , FIO. 5 , lOX. ' STD  -',F10  5/lOX 
*'MIN  OF  ET  -',F10.5,10X,'MAX  OF  ET  -', FIO. 5/lOX, 'LOWER  LIMIT  =' 
*F10.5,10X,'UPPER  LIMIT  -' , FIO. 5/lOX, '#  OF  2STD  CONTAIN  TRUE  VALUE 
*-', FIO. 5, lOX,' PERCENT  OF  2STD  CONTAIN  TRUE  VALUE  -'  FIO  5/lOX 
*'#  OF  3STD  CONTAIN  TRUE  VALUE  -', FIO. 5 , lOX, ' PERCENT  OF  3STD  CONTAI 
*N  TRUE  VALUE  -', FIO. 5/lOX, 'MEAN  OF  STD  OF  EXPECTED  TIME  TO  MATURIT 
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*Y  -',F10.5/10X,' VARIANCE  OF  STD  OF  EXPECTED  TIME  TO  MATURITY  =' , 
*F10.5,10X,'STD  OF  STD  OF  EXPECTED  TIME  TO  MATURITY  ='  FIG  5) 
WRITE(6,116) 

116  FORMATC-' ,5X,'THE  FOLLOWING  VALUE  IS  WHEN  THE  PROBABILITY  IS  ZER 
*0:  ') 

WRITE (6, 117)  TTO,  ETO(KK) ,  SUMO,XOBAR,ETSQO,VARO, STDO,ETOMIN, 
*  ETOMAX , DLOWO , UPPERO , C0UN02 , PER02 , COUN03 , PER03 , STOBAR , VAROST 
*,STDOST 

117  FORMAT('-M0X,'TT=' ,  FIG.  5,  lOX,' EXPECTATION  OF  TIME  TO  THE  MATURI 
*TY  =',  F1G.5/1GX,'SUM-',F1G.5,10X,'XBAR-',F10.5  /IGX/SUM  OF  SQUAR 
*E  =', FIG. 5, /lOX,' VARIANCE  -' ,F10. 5 , lOX, ' STD  =' ,F10 . 5/lGX, 

*'MIN  OF  ET  -' ,F10.5,10X,'MAX  OF  ET  -', FIG. 5/lOX, 'LOWER  LIMIT  =' 
*F1G.5,1GX,'UPPER  LIMIT  =' , FIG. 5/lOX, '#  OF  2STD  CONTAIN  TRUE  VALUE 
*=  ', FIG. 5, lOX,' PERCENT  OF  2STD  CONTAIN  TRUE  VALUE  =',F10.5/1GX 
*,'#  OF  3STD  CONTAIN  TRUE  VALUE  =', FIG . 5 , IGX, ' PERCENT  OF  3STD  CONTA 
*IN  TRUE  VALUE  =', FIG. 5/lGX, 'MEAN  OF  STD  OF  EXPCTED  TIME  TO  MATURI 
*TY  =', FIG. 5/lOX,' VARIANCE  OF  STD  OF  EXPECTED  TIME  TO  MATURITY  =' 
*,F1G.5,1GX,'STD  OF  STD  OF  EXPECTED  TIME  TO  MATURITY  ='  FIG  5) 
6   CONTINUE 

WKVEC(1)=1.DG 
GO  TO  5G1 
5GG   STOP 
END 
C 

c         ************************************ 
c  *  * 

C         *        DEFINE  FN  FUNCTION        * 
C        *  * 

c  ************************************ 

c 

FUNCTION  FN(X,Y) 

IMPLICIT  REAL*8  (A-H,0-Z) 

COMMON   UX,  UY,  VX,  VY,  PXY 

PP-PXY*PXY 

A-DSQRT(l-PP) 

B-(X-UX)/VX 

C-(Y-UY)/VY 

D  -  (B*B)  -  (2*PXY*B*C)  +  (C*C) 

FN  -  X  *  (1/(2*3. 14159*VX*VY*A))*DEXP((-1/(2*(1-PP)))*D) 
RETURN 

END 

C 

c  **************************** 

c  *  * 

C  *     DEFINE  FD  FUNCTION   * 

C  *  * 

c  **************************** 

c 
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FUNCTION  FD(X,Y) 

IMPLICIT  REAL*8  (A-H,0-Z) 

COMMON  UX,  UY,  VX,  VY,  PXY 

PP=PXY*PXY 

A=DSQRT(1-PP) 

B=(X-UX)/VX 

C=(Y-UY)/VY 

D  =  (B*B)  -  (2*PXY*B*C)  +  (C*C) 

FD  =  (1/(2*3. 14159*VX*VY*A))*DEXP((-1/(2*(1-PP)))*D) 

RETURN 

END 
C 

c        ************************************ 
c       *  * 

C         *      DEFINE  LOWER  LIMIT  * 

C         *  ON  Y  * 

C        *  * 

c        ************************************ 

c 

FUNCTION  AY(X) 

IMPLICIT  REAL*8  (A-H,0-Z) 

AY-X 

RETURN 

END 
C 

c      ******************************** 
c     *  * 

C      *     DEFINE  UPPER  LIMITS      * 
C      *         ON  Y  * 

C      *  * 

c      ******************************** 

c 

FUNCTION  BY(X) 
IMPLICIT  REAL*8   (A-H,0-Z) 
BY-30.D0 
RETURN 
END 
/* 

//LKED.SYSLIB   DD 

//  DD 

//  DD 

//  DD  DSN-SYSl.IMSL.SPFLIB.DISP-SHR 

//GO. SYS IN   DD  * 

6  9  2  2  0 

6  9  2  2  25 

6  9  2  2  50 

6  9  2  3  0 

6  9  2  3  25 
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TABLE  1 
simulate  100  times 


/i  =  6,   a  =2 


P   =  0.0 


equally  spaced     unequally  spaced 


^y       °y  "^  -  15  .     30         15      30 


true  value  5.63 

A 

estimated  T         5.74      5.70      5.77    5.73 
confidence  (5.50,     (5.51,     (5.53,   (5.54, 

A 

interval  of  T        5.98)     5.90)     6.01)    5.92) 

A 

mean  of  std(T)      1.59      1.16      1.53    1.11 
true  value  5 .  62 

A 

estimated  T         5.37      5.62      5.38    5.63 
confidence  (5.06,     (5.42,     (5.09,   (5.44 

A 

interval  of  T        5.68)     5.81)     5.68)    5.82) 

A 

mean  of  std(T)      1.68      1.19      1.62    1.15 
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TABLE  1  (continued) 
simulated  100  times 


^  =  6,     a  =  2 
'^x         X 


p   =  0.0 

equally  spaced     unequally  spaced 


My   ffy  n  -   15       30       15       30 


true  value  5.89 

A 

estimated  T         5.94      5.95     5.94     5.94 

3    

confidence  (5.83,     (5.89,    (5.83,    (5.88, 

A 

interval  of  T        6.05)     6.02)    6.05)    6.01) 

A 

mean  of  std(T)      0.49      0.34     0.53     0.37 
12  

true  value  5 . 85 

A 

estimated  T        6.06      6.05     6.06     6.05 

4    - 

confidence  (5.97,     (5.98,    (5.97,    (5.98, 

A 

interval  of  T        6.15)     6.11)    6.16)    6.13) 

A 

mean  of  std(T)      0.45      0.34     0.48     0.36 
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TABLE  1  (continued) 
simulated  100  times 


^x    '    X 


P   =  0.0 

equally  spaced     unequally  spaced 


/^y   «^y  n  =    15       30       15       30 


true  value  5.96 


15 


estimated  T  6.02      6.03      6.04    6.05 

confidence  (5.94,     (5.97,     (5.96,   (5.99, 

A 

interval  of  T  6.09)     6.08)     6.12)    6.11) 

A 

mean  of  std(T)  0.38      0.27      0.43    0.31 

true  value  5.93 

estimated  T  6.02      6.03      6.01    6.01 

confidence  (5.95,     (5.97,     (5.92,   (5.95, 

A 

interval  of  T  6.10)     6.08)     6.09)    6.06) 

A 

mean  of  std(T)  0.37      0.27      0.42    0.30 
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TABLE  1  (continued) 
simulated  100  times 
M  =  6,    a  =  2 

X  X 


p   =  0.25 


equally  spaced     unequally  spaced 


^y        "y  n  =  15   -     30         15      30 


true  value  5.75 


estimated  T         5.54      5.65      5.56    5.69 
confidence  (5.29,     (5.48,     (5.32,   (5.53, 

A 

interval  of  T        5.79)     5.82)     5.79)    5.85) 

A 

mean  of  std(T)      1.65      1.20      1.59    1.16 
true  value  5.77 

A 

estimated  T         5.55      5.61      5.55    5.64 
confidence  (5.30,     (5.41,     (5.30,   (5.45, 

A 

interval  of  T        5.80)     5.81)     5.80)    5.82) 

A 

mean  of  std(T)      1.58      1.17      1.52    1.12 
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TABLE  1  (continued) 
simulated  100  times 


M  =  6,   a     =  2 

X      '      X 


p  =  0.25 


equally  spaced     unequally  spaced 


/^y   '^y  n  -      15  .     30        15      30 


true  value  5.96 


12 


estimated  T         5.98      5.97      5.99    5.97 
confidence  (5.87,     (5.90,     (5.89,   (5.89, 

A 

interval  of  T        6.08)     6.04)     6.09)    6.04) 

A 

mean  of  std(T)      0.50      0.35      0.53    0.37 
true  value  5  .  94 

A 

estimated  T         5.90      5.93      5.91    5.94 
confidence  (5.80,     (5.87,     (5.80,   (5.87, 

A 

interval  of  T        6.00)     6.00)     6.01)    6.01) 

A 

mean  of  std(T)      0.53      0.36      0.55     0.38 
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TABLE  1  (continued) 
simulated  100  times 


M  =  6,    a  -  2 


P  =  0.25 

equally  spaced     unequally  spaced 


/*y   «^y  n  =     15  .    30        15      30 


true  value  5.99 


15 


estimated  T         5.94      5.94      5.92    5.95 
confidence  (5.86,     (5.88,     (5.84,   (5.89, 

A 

interval  of  T        6.02)     5.99)     6.01)   6.00) 

A 

mean  of  std(T)      0.37      0.26      0,42     0.30 
true  value  5.98 

A 

estimated  T         6.01      6.01      6.00    6.03 
confidence  (5.93,     (5.96,     (5.91,   (5.97, 

A 

interval  of  T       6.09)     6.06)     6.08)   6.09) 

A 

mean  of  std(T)      0.38      0.27      0.43    0.30 
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TABLE  1  (  continued  ) 

simulated  100  times 

M  =  6,   a  -  2 
^x     '    X 

P   =  0.5 

equally  spaced     unequally  spaced 


A'y   t^y  n  =     15  ■     30         15      30 


true  value  5.87 

A 

estimated  T        5.81      5.67      5.79    5.66 

2    - 

confidence         (5.57,     (5.49,     (5.55,   (5.48, 

A 

interval  of  T        6.06)     5.86)     6.03)    5.83) 

A 

mean  of  std(T)      1.46      1.15      1.42    1.10 
9 

true  value  5.92 

A 

estimated  T         5.58      5.61      5.56    5.60 

3    

confidence  (5.32,     (5.43,     (5.30,   (5.42, 

A 

interval  of  T        5.84)     5.80)     5.82)    5.78) 

A 

mean  of  std(T)      1.55      1.17      1.50    1.12 
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TABLE  1  (continued) 
simulated  100  times 


M  -  6,    a  =  2 

X  X 


P   =  0.5 

equally  spaced     unequally  spaced 


^y   '^y  n  =   15   -     30         15      30 


true  value  6.00 

A 

estimated  T        5.82      5.94      5.79    5.94 

3    - 

confidence  (5.71,     (5.88,     (5.67,   (5.87, 

A 

interval  of  T        5.93)     6.01)     5.91)    6.01) 

A 

mean  of  std(T)      0.54      0.36      0.57    0.38 
12  

true  value  6.01 

A 

estimated  T         5.97      5.98      5.96    5.95 

4 

confidence  (5.88,     (5.92,     (5.87,   (5.87, 

A 

interval  of  T        6.07)     6.05)     6.06)    6.02) 

A 

mean  of  std(T)      0.46      0.33      0.50    0.35 
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TABLE  1  (continued) 
simulated  100  times 


M  =  6,    a     =   1 

X  X 


p   =  0.5 

equally  spaced     unequally  spaced 


l^y       °y  n  -    15  .     30        15      30 


true  value  6.01 


15 


estimated  T  5.99      6.00      5.96    5.98 

confidence  (5.92,     (5.95,     (5.88,   (5.93, 

A 

interval  of  T  6.07)     6.05)     6.05)   6.03) 

A 

mean  of  std(T)  0.37      0.27      0.42    0.30 

true  value  6 . 02 

A 

estimated  T  5.96      6.01      5.95    5.98 

confidence  (5.89,     (5.95,     (5.87,   (5.91, 

A 

interval  of  T  6.03)     6.07)     6.04)    6.05) 

mean  of  std(T)  0.37       0.27     0.42     0.30 


simulated  100  times 
M^  =  6,    a^  =  2 
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TABLE  2 


A*     cr 

y   y 

0.0 

0.25 

0.5 

true  value 

5.63 

5.75 

5.87 

n  -  15 

(0.84) 

(0.87) 

(0.83) 

equally 

[0.91] 

[0.95] 

[0.92] 

spaced 

30 

(0.89) 

(0.95) 

(0.93) 

[0.91] 

[0.98] 

[0.96] 

2    

15 

(0.84) 

(0.91) 

(0.85) 

unequally 

[0.94] 

[0.96] 

[0.93] 

spaced 

30 

(0.89) 

(0.96) 

(0.94) 

9   - 

[0.92] 

[0.97] 

[0.96] 

true  value 

5.62 

,5.77 

5.92 

15 

(0.88) 

(0.86) 

(0.93) 

equally 

[0.94] 

[0.97] 

[0.99] 

spaced 

30 

(0.99) 

(0.95) 

(0.96) 

3    

[0.99] 

[0.96] 

[0.97] 

15 

(0.91) 

(0.90) 

(0.97) 

unequally 

[0.96] 

[0.98] 

[0.99] 

spaced 

30 

(0.98) 

(0.95) 

(0.96) 

[0.99] 

[0.96] 

[0.99] 
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TABLE  2  (continued) 


simulated  100  times 


^x  =  6-    -X  =  2 


u         a 

y      y 

0.0 

0.25 

0.5 

true  value 

5.89 

5.96 

6.00 

n  -  15 

(0.89) 

(0.91) 

(0.97) 

equally 

[0.99] 

[0.99] 

[1.00] 

spaced 

30 

(0.94) 

(0.94) 

(0.98) 

3 

[1.00] 

[0.99] 

[1.00] 

15 

(0.92) 

(0.93) 

(0.96) 

unequally 

[1.00] 

[0.99] 

[1.00] 

spaced 

30 

(0.95) 

(0.93) 

(0.97) 

L2 

[1.00] 

[0.97] 

[1.00] 

true  value 

5.85 

5.94 

6.01 

15 

(0.90) 

(0.91) 

(0.96) 

equally 

[0.98] 

[0.99] 

[1.00] 

spaced 

30 

(0.86) 

(0.96) 

(0.95) 

[0.97] 

[1.00] 

[1.00] 

4 

15 

(0.89) 

(0.94) 

(0.93) 

unequally 

[0.97] 

[1.00] 

[0.98] 

spaced 

30 

(0.85) 

(0.97) 

(0.95) 

[0.97] 

[1.00] 

[0.99] 

simulated  100  times 


M^=6,    a^=2 
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TABLE  2  (continued) 


^i  a 

y       y 


0.0 


0.25 


unequally 
spaced 


15 


30 


(0.91) 
[0.99] 
(0.97) 
[0.99] 


(0.94) 
[0.99] 
(0.96) 
[1.00] 


0.5 


true  value 

5.96 

5.99 

6.01 

n  -  15 

(0.95) 

(0.93) 

(0.92) 

equally 

[1.00] 

[1.00] 

[1.00] 

spaced 

30 

(0.93) 

(0.94) 

(0.97) 

[0.98] 

[0.99] 

[1.00] 

4   

15 

(0.93) 

(0.93) 

(0.94) 

unequally 

[1.00] 

[0.98] 

[0.99] 

spaced 

30 

(0.93) 

(0.96) 

(0.97) 

15  

[0.99] 

[0.99] 

[1.00] 

true  value 

5.93 

5.98 

6.02 

15 

(0.91) 

(0.95) 

(0.97) 

equally 

[1.00] 

[1.00] 

[1.00] 

spaced 

30 

(0.94) 

(0.97) 

(0.92) 

[1.00] 

[1.00] 

[1.00] 

5    

(0.92) 
[1.00] 
(0.90) 
[1.00] 
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The  value  in  the  parenthesis  [bracket]  is  the  percent  of  2  std  C.I.'s 
[3  std  C.I.'s]  which  contain  the  true  value. 
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ABSTRACT 
Consider  an  insect  that  progresses  through  developmental  stages 
throughout  its  lifetime.  We  here  consider  only  one  transition,  i.e., 
from  non-maturity  to  maturity.  We  take  independent  samples  of  fixed 
size  and  record  the  numbers  of  both  mature  and  non-mature  organisms, 
noting  mortality  in  each  group,  at  various  observing  times.  We 
estimate  the  average  time  to  maturity,  and  its  corresponding  standard 
error.  We  also  present  a  simulation  study  comparing  the  estimated  time 
to  maturity  with  its  theoretical  value. 


